%函数：get_result(初始概率*100，概率开始增加的抽数，每次增加的概率*100，保底抽数)
%输出值：返回抽中的期望抽数和每一抽的期望概率
%此函数只适用于前N次概率不变，第K次开始概率均匀增加的抽卡模型，且第M次时概率直接变为1

function [num,pro] = get_result2(p_first,inflexion,p_add,ends)
%初始概率，%概率转折点，%每次增加概率，%保底抽数
p_add_num=ends-inflexion;%概率变化的抽数
A=ones(1,inflexion-1)*p_first/100;
B=linspace(p_add/100,p_add_num*p_add/100,p_add_num)+p_first/100;
P=[A,B,1];%第K次抽中的概率
Pi=ones(1,ends);a=1;q=0;
for k=1:ends
    Pi(k)=a*P(k);%恰好第K次抽中概率
    q=q+Pi(k)*k;%出货期望抽数
    a=a*(1-P(k));
end
plot(Pi);
format long
num = q;
pro = 1/q*100;
format default